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Abstract 

Based on ab initio calculation, we propose a new 
structure for the fundamental excitation of the re- 
constructed 30° partial dislocation in silicon. This 
soliton has a rare structure involving a five-fold 
coordinated atom near the dislocation core. The 
unique electronic structure of this defect is con- 
sistent with the electron spin resonance signature 
of the hitherto enigmatic thermally stable R cen- 
ter of plastically deformed silicon. This identifica- 
tion suggests the possibility of an experimental de- 
termination of the density of solitons, a key defect 
in understanding the plastic flow of the material. 
PACS 71.15.Mb 

Dislocations are central to the understanding of the 
mechanical response of materials. The mechanical 
behavior of any crystalline material is determined by 
a hierarchy of crystalline defects of successively lower 
dimension. Grain boundaries are two dimensional de- 
fects that control the evolution of the microstructure 
of the material. The creation and motion of disloca- 
tions, which are one dimensional extended topological 
defects of the lattice, mediate the plastic response of 
a crystal to external stress. In silicon, which has a 
bipartite lattice, the primary mobile dislocations are 
the screw and the 60° dislocation, which dissociate 
into more primitive one-dimensional partial disloca- 
tions bounding stacking faults. The mobility of dislo- 
cations in high Peierls barrier materials such as silicon 
is effected by the motion of complex zero-dimensional 
local defects known as kinks, where the dislocation 
center skips from one row of atoms to a neighbour- 
ing row. The low energy kinks along the 30° partial 
dislocation in silicon have been shown to also involve 
a composite structure where kinks bind with zero- 
dimensional soliton excitations in the reconstructed 
ground state of the dislocation core. These soliton 
excitations are also known as "anti-phase defects" 



(APDs) 0. 

Here we report the results of an ab initio study ex- 
ploring the lattice and electronic structures, excita- 
tion energy, and the density of these APDs which are 
the simplest, lowest energy, fundamental excitations 
of the dislocations in the hierarchy ultimately leading 
to the macroscopic mechanical behavior of the solid. 
The soliton is associated with an atom in the disloca- 
tion core which is not part of a reconstructed dimer. 
In the simple, conventional picture, this atom (hence- 
forth to be referred to as the "soliton atom" ) only has 
three bonds and therefore an unpaired electron. This 
simple model, however, does not lead to predictions 
consistent with any of the observed electron spin res- 
onance (ESR) signals associated with plastically de- 
formed silicon. 

We propose a new structure of the soliton. We 
find that the ground state of the soliton has an un- 
expected structure with electronic states which are 
consistent with the most stable ESR center in plas- 
tically deformed silicon, the only one which remains 
after careful annealing. The reason why the natu- 
ral connection between this center and the lowest en- 
ergy excitation of the dislocation core has not been 
made previously is that the observed ESR center has 
a highly unusual symmetry. In support of our theory 
for the structure of the soliton, we gather here several 
pieces of evidence from both reports of ESR results 
and our own ab initio calculations. The final com- 
bined ab miiio-experimental identification which we 
make allows for the possibility that future more pre- 
cise measurements of the ESR center density could be 
used to determine experimentally the soliton density, 
a key physical quantity in the process of deformation. 

Figure |^ reviews the basic geometry of the 30° 
partial dislocation studied in this work. The dislo- 
cation is the one-dimensional boundary defining the 
edge of a half- planar (111) stacking fault. Atoms in 
the central core of the dislocation (shaded gray in 
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Figure 1: (a) Three dimensional view of a 30° partial 
dislocation, with a soliton domain wall, (b) schematic 
drawing of the same in the (111) plane. The grey 
circles represent the atoms in the central core of the 
dislocation. 



the figure) are connected to the bulk with only three 
bonds per atom. The dislocation undergoes a recon- 
struction whereby the core atoms pair up in dimers 
forming intra-core bonds and thus become four-fold 
coordinated. This reconstruction breaks the original 
translational symmetry along the core and doubles 
the primitive repeat distance along the core axis. As- 
sociated with this broken translational symmetry is 
a low energy soliton defect where, by creating a sin- 
gle unpaired core atom, the system may change along 
the line from one of the two symmetry related degen- 
erate ground state phases to the other. Because the 
soliton atom is expected to have a dangling bond, 
it is natural to look for an ESR signal for this de- 
fect. The relatively low energy we expect for such an 
excitation leads us to expect a relatively large equi- 
librium population at silicon annealing temperatures 
(~ 900K) and therefore that the ESR signal would 
not anneal out as quickly as other signals associated 
with the formation and motion of dislocations. We 
further would expect this signal to be detected in all 
systems which contain 30° partial dislocations. 

Indeed, it was discovered over thirty years ago that 
plastically deformed silicon gives a wide variety of 
ESR signals j^. (For an excellent review, see [^.) 
Out of the dozens of ESR centers, four have been 
identified as associated with the 30° dislocation core. 
The literature refers to these centers as Si-Kl, Si-K2, 
Si-Y and Si-R. The Kl and K2 defects have been 
identified to be electronic excitations of the same 
structural defect. Kisielowski-Kemmerich made 
the currently accepted identification of the K and Y 
defects. It is well known that the first three of the 



aforementioned ESR signals anneal out over the (tem- 
perature dependent) time-scale of about an hour |3]. 
Only one signal remains, the one labelled R |, |]. 
This center is "thermally stable" (it does not anneal 
out) and is observed even at high deformation tem- 
peratures (> 900° K) where the other signals anneal 
out too quickly to be observed. The R signal is the 
residue of Y after annealing, and is very similar to Y 
in its diminished anisotropy and large width js] . 

Many authors have studied excitations of disloca- 
tions in silicon |, ||. The nature of the ESR 
centers and the soliton excitation energy of the 30° 
partial dislocation, however, have yet to be addressed 
with modern ab initio techniques. In order to investi- 
gate the electronic structure of low energy excitations 
of the 30° partial dislocation core, we embarked upon 
a density functional study of the system. 

To prepare approximately relaxed initial ionic con- 
figurations with the correct bonding topology, we first 
relaxed lattices containing dislocation cores using 
the Stillinger- Weber (SW) inter- atomic potential [p^. 
While doing this, we discovered that the soliton atom 
moves out of line with respect to the dislocation core. 
To probe this interesting possibility further, we car- 
ried out calculations within the plane wave total en- 
ergy density functional approach [pT|. To describe the 
electron-electron interactions we used the Perdew- 
ZungerfT^ parameterization of the Ceperly- Alder [T^ 
exchange-correlation energy of the uniform electron 
gas. To describe the electron-ion interactions we 
used a non-local pseudopotential of the Kleinmann- 
Bylander form|Q. The electronic wave functions 
were expanded in a plane- wave basis up to a cutoff of 
8Ry. 

All super-cells used in this study have the same size 
in the plane perpendicular to the (110) dislocation 
axis. Two partial dislocations of equal but opposite 
Burgers vectors at a separation of 14 A cut through 
this plane. Following Bigger et al. |jl^, the lattice 
vectors are arranged so that the periodic dislocation 
array has a quadrupolar arrangement. Each cell con- 
tains forty-eight atoms per bilayer stacked along the 
dislocation core direction. To calculate the excitation 
energy of the soliton, it is also necessary to calculate 
the energy of the perfectly reconstructed dislocation. 
However, the smallest super-cell which is commensu- 
rate with both structures contains six bilayers (288 
atoms). It is possible, however, to reduce the compu- 
tational time by using two different super-cells. For 
the reconstructed case, the super-cell contains two 
bilayers along the dislocation line, while the soliton 
structure contains three, (96 and 144 atoms, respec- 
tively) . The lattice vectors were obtained by relaxing 
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Figure 2: Coplanar atoms near the soliton: (a) in 
the proposed ground state for the sohton, (b) con- 
ventional structure with the sohton atom in hne with 
the remaining atoms in the dislocation core. The no- 
tation for the atoms is the same as in Fig. |l|: atom 
A is the soliton atom, B is the five-fold coordinated 
bulk atom (see text), C is a bulk atom bonded to B 
opposite from A. 

a completely reconstructed dislocation within the SW 
model in the 96 atom cell. The three bilayer cell was 
then obtained from this by multiplying the lattice 
vector that points along the dislocation axis by |. 

To minimize errors from k-point sampling, ba- 
sis set truncation and super-cell effects, we com- 
pute differences of energy differences, as follows. For 
each super-cell, we generated a completely unrecon- 
structed configuration where all the core atoms have 
only three bonds. These structures can be realized 
in both super-cells, and thus serve as the reference 
point. The final excitation energy is the difference 
between the deviations in the energy from the unre- 
constructed structure in each cell. By keeping the 
lattice vectors fixed throughout the calculations, we 
more closely simulate the environment which widely 
separated solitons would experience along a recon- 
structed dislocation line. To test for the impact of 
strain effects, we repeated all calculations using the 
corresponding lattice vectors of a bulk silicon system 
at the ab initio lattice constant. Our results did not 
change significantly. 

To ensure maximum transferability of results be- 
tween the two cells, the calculations employed k- 
point sets which give identical sampling of the Bril- 
louin zone for the two super-cells: {(0,0, ±1/4)} for 
the three bilayer ceU and {(0, 0, ±1/6), (0, 0,1/2)} for 
the two bilayer cell. To find relaxed structures, we 
moved all ions along the Hellmann-Feynman forces 
until the ionic forces were less than 0.02 eV/A. Typ- 
ically, this was accomplished in 40 ionic steps, where 
between ionic steps we made 10-15 electronic re- 
laxation steps using the analytically continued func- 
tional approach [p^. 

Figure ^ shows the projection in the (110) plane 



of our ab initio results for the structure of the soli- 
ton. (The same atoms arc shown in Figure |l| in full 
structural context.) The structure on the left is our 
prediction for the ground state. In this configura- 
tion, there is a five-fold coordinated bulk atom (B) 
in the immediate neighbouring row to the disloca- 
tion core. Its new, fifth neighbor is the soliton atom 
(A) . Our ab initio results show that the conventional 
structure on the right (generated by keeping the soli- 
ton atom coUinear with the dislocation core) is not 
only higher in energy but also spontaneously decays 
into the ground state on the left. 

The ah initio excitation energy of the soliton is 
0.65 ± (« 0.2) eV, where we attribute most of 
the uncertainty to the uncontrolled local density ap- 
proximation and super-cell effects. This energy corre- 
sponds to a density p = ^e~^/^'^ in the range of 10~^ 
to 2 X 10"'^ solitons per core atom at 900 K, which 
is consistent with the observed densities of the ESR 
centers: the R center represents a "fraction" of the 
density of the Y center, which is estimated in the ex- 
perimental literature to be about 0.01 per core atom 
(Table II in §, Table 2 in 0, and §). (The factor 
of one half in p comes from the fact that in a given 
phase of the ground state, only half of the core atoms 
represent possible soliton sites.) Given the exponen- 
tial sensitivity of the density, we find this agreement 
encouraging, particularly as the energy of soliton is 
very low compared typical point defects in silicon. 

A great advantage of ab initio calculations, beyond 
their reliability, is that they also yield the electronic 
states, in particular giving information about their 
spatial symmetry. In Figure ^ wc plot the angu- 
lar momentum decomposition of the local density of 
states as obtained from the Kleinmann-Bylander pro- 
jections of the electronic eigenstates. Panel (a) shows, 
for an atom far from the core, the familiar concentra- 
tion of p-like states at the top of the valence band. 
To explore the nature of the soliton state, we compare 
this to the local densities of states for the soliton atom 
in the proposed (^) and conventional (^) configura- 
tions. We also plot the local density of states for the 
quasi-fivefold coordinated atom (|^d). The appear- 
ance of the peak near the top of the valence band in 
the s channel of the soliton atom (and the correspond- 
ing diminution in the p channel) in its ground state 
) shows that the state associated with the soli- 
ton is much less anisotropic than the simple dangling 
p-like bond on the soliton atom in the conventional 
picture (|^c). We further note an enhancement at the 
same energy in the s channel of the quasi-fivefold co- 
ordinated atom (^), which indicates that the un- 
paired electron is shared between this atom and the 
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Figure 3: Local density of states, calculated by act- 
ing on the filled bands with the Kleinmann-Bylander 
projectors centered at (a) an atom deep in the bulk 
, (b) the soliton atom in its proposed ground state, 

(c) a conventional soliton atom with a dangling bond, 

(d) the quasi five-fold coordinated atom. The hori- 
zontal axis is the energy (eV), the scale of the vertical 
axis is arbitrary but the same for all four plots. Solid 
and dashed lines represent densities in the s and p 
channels respectively. 



soliton atom. Defects in the dislocation core there- 
fore need not be associated with strongly directional 
electronic states, as has been previously assumed in 
identifications of ESR centers. 

There is a direct connection between the symme- 
try of the electronic state of an unpaired spin and 
the symmetry of the corresponding ESR signal as de- 
scribed by its effective g tensor. The off diagonal 
elements of this tensor involve a sum over matrix el- 
ements of the form ((/)o|Li|0e)(0e|Lj |0o)j where 4>q is 
the unpaired state and (pe are the excited states [p^. 
In general 0o can be broken into angular momen- 
tum components (as in Figure |3|), of which the s 
wave component makes no contribution to the pre- 
ceding matrix elements. In general, then, we expect 
the anisotropy of g to be proportional to the popu- 
lation of the p channel. (Higher angular momentum 
components are negligible for filled states in silicon.) 
This population, the area under the peak associated 
with the unpaired electron in the p channel, drops by 
about a factor of two as the soliton moves from its 
symmetrical dangling bond configuration (^) to our 
proposed state (||b). The literature contains qualita- 



Figure 4: Two dimensional slice of the total charge 
density from our ab initio calculation, through the 
plane containing the atoms A, B, C, as labelled in 
the previous figures. The new A-B bond is nearly as 
pronounced as the bulk B-C bond, which is clearly 
weakened compared to the other (vertical) bulk bond 
of atom B seen in the figure. 



five observations of the decreased anisotropy of the R 
signal and one quantitative comparison which comes 
from measurements of the Y signal, of which the R 
is the residual after annealing. In [Q, this anisotropy 
is compared directly with that of the Kl,2 centers, 
which have typical dangling bonds, and is shown to 
be less by about a factor of two, in agreement with 
our electronic structure results. 

To explore the nature of the bonding near the soli- 
ton atom, we plot the total valence charge density 
in Figure |^, which shows that the soliton atom (A) 
makes a weak bond with the neighbouring five-fold 
coordinated atom in the bulk (B). The new bond (A- 
B) of the bulk atom is very similar to a now weak- 
ened but previously existing bond in the spatially 
opposite direction (B-C). Such five-fold coordinated 
structures have been considered previously in silicon 
by Pantelides (lO) and more recently by Duesbery et 
al. . Through ab initio studies it was demonstrated 
that such five- fold defects in amorphous silicon should 
show similar anomalies in the angular momentum de- 
composition of the local density of states ||2l]|. 

In conclusion, we have presented ab initio results 
indicating that the ground state of the soliton has 
an unusual structure involving a five-fold coordinated 
atom and a correspondingly unusual electronic struc- 
ture. The excitation energy we calculate for this new 
defect corresponds to a thermal equilibrium density 
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which is compatible with the observed density of the 
R center, which is the only thermally stable param- 
agnetic center associated with the 30° partial dislo- 
cation. In line with our notion of the soliton being 
the fundamental excitation of the reconstructed dis- 
location core, the R center is observed independent of 
the method of deformation and in proportion to the 
dislocation density. Our calculations show that the 
electron state of the soliton has a reduced anisotropy 
compared to that of a simple dangling bond which 
corresponds in magnitude to the puzzling reduction 
in anisotropy of the ESR signal of the R center. In the 
groimd state structure which we propose, the soliton 
atom makes a weak bond with a neighbouring bulk 
atom and thus gives rise to an amorphous-like bond- 
ing arrangement. This could explain in part the oft 
noted similarity of the ESR signature of the R cen- 
ter to that of amorphous silicon. Based on the above 
arguments and results, we propose that the domain 
walls in the reconstruction of the 30° partial disloca- 
tion and the R centers observed in ESR experiments 
are one and the same. Any viable competing theory 
which does not identify the R signal with the soliton 
must both predict a more plausible intrinsic excita- 
tion corresponding to the the R center and explain 
why the unpaired electron of the low energy soliton 
is not detected in ESR experiments. 
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